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I. INTRODUCTION 



(N 
> 

, By studying transitions of higher dimensional nonintegrable maps to integrable ones, we showed in our previous 

o. fl 

paper 1] that the periodicity conditions, which are 'uncorrelated' in generic systems, become 'fully correlated' in the 
I/"") , integrable maps and one gets full manifolds, instead of isolated points, characterized by invariants of the map. The 



purpose of this paper is to establish a criterion which distinguishes the fully correlated and uncorrelated periodicity 



conditions. We will also discuss in detail how the transition takes place between them. 



We consider a rational map on C d 

x=(xx,x 2 ,...,x d ) - X=(X 1 ,X 2 ,...,X d )=xM. (1) 



1 We are interested in the behaviour of the sequence: x — * X^ — > X^ In the study of the sequence of the 

map, the periodic points of the map decide the most important nature of the sequence. To sec their role let us begin 
our discussion by introducing some words. The precise meaning of them will be explained in §3. 

Let us assume that the map has p invariants. The periodicity conditions are called 'uncorrelated' if they determine 
all positions of the periodic points dependent on the values of the invariants. They are called 'fully correlated' if they 
do not determine the position of the periodic points but impose relations on the values of the invariants. An 'invariant 
variety of periodic points' is a variety determined by the invariants of the map alone, such that every point on the 
variety can be an initial point of the iteration, which stays on the variety before it returns to the same point. We 
have shown in [l| the following propositions: 
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• A set of periodic points form a set of isolated points if the periodicity conditions are uncorrelated, and an 
invariant variety of dimension p if they are fully correlated. 

• An invariant variety of periodic points and a set of isolated periodic points do not exist in one map simultaneously. 

We would like to prove the following theorem in this paper, from which the above second proposition immediately 
follows. 

Theorem: 

• When p < d/2, there is no set of fully correlated periodicity conditions. 

• When p > d/2, sets of correlated and uncorrelated periodicity conditions do not exist in one map simultaneously. 

We don't know exactly the correspondence between the existence of an invariant variety of periodic points and 
integrability of the map. Nevertheless we have presented in [J many examples of invariant variety found in some 
known integrable rational maps. On the other hand, the isolated periodic points, which characterize non-integrable 
maps, are derived by solving uncorrelated periodicity conditions. Therefore it will be worthwhile to study the 
transition of periodicity conditions from uncorrelated regime to fully correlated one. This will be our second topic 
in this paper. We would like to discuss in detail a simple example and see how the transition takes place when a 
parameter introduced to the map approaches from one regime to the other. In particular we focuss our attention to 
the behaviour of the Julia set, the closure of the set of repulsive periodic points of a map. 

This paper will be organized as follows. We study in §2 a general scheme of higher dimensional maps with p 
invariants. The scheme enables us to investigate integrable and nonintegrable maps on an equal footing. The proof of 
the theorem will be given in §3. Based on the results of the previous sections we develop in §4 a method generating 
invariant varieties of periodic points for some simple cases. We study in §5 how the transition takes place between 
the uncorrelated and fully correlated periodicity conditions. 

II. DIMENSIONAL REDUCTION BY INVARIANTS 

The integrability of a dynamical system with continuous time evolution is guaranteed by the existence of sufficient 
number of invariants of motion due to the Liouville- Arnold theorem. There is, however, no corresponding theorem in 
the case of iteration of maps. In this section we study the iteration of a map Q on C d with p invariants. 

To begin with we study a general scheme of a map with p invariants Hi(x), i — 1,2,..., p. If the functions 
fj( x )i j = 1,2,..., d satisfy 

H i (fi(x),---,f d (x))=H i (x 1 ,x 2 ,---,Xd),, * = l,2,...,p (2) 
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the d dimensional map 

(x 1 ,x 2 ,---,Xd) -» (X 1 ,X 2 ,---,X d ) = (f 1 (x),f 2 (x),---,fd(x)\ (3) 

restricts to any common level set of H±, ■ ■ ■ , H p . More generally, the invariants Hi may be specified implicitly, for 
example by 

Vi(h{x),h(x),->>,U{x),H?)=Q, i=l,2,...,p. (4) 

In general the equations Q fix a level set V(h), 

V(h) = {z| Hi{x) =hi, i= 1, 2, (5) 

with h\, h 2 , hp the values of the invariants determined by the initial data of the map. V(h) is a d — p dimensional 
variety in C d and © defines a map whose trajectory of iteration is constrained to the variety. Different choices of f\ 
correspond to different parametrizations of the trajectory. 

As an example, let H(x, y) — x + y, i.e, 

fix, y,H) = x + y - H. (6) 

If we fix fi(x,y) by 

fi(x,y)=xy, (7) 

then, by solving fi + f% = x + y, we find f 2 (x, y) = x + y — xy, and the map 

(x,y) -» (X,Y) = (xy, x + y-xy^j (8) 

restricts to any line x + y = h. Using this fact we can reduce the map (JSJ) to one dimension, 

X = x(h-x), Y = h-X. 

We see that the 2-dim map (JSJ is equivalent to a logistic map, hence is not integrable although it has the invariant 
H = x + y. 

Next, consider an example where the level set is a hyperbola: 

<p(x, y, H) = (1 - bx)y - (1 - cx)H, (9) 

and fix the same function (7J for fi(x, y). We obtain 

f 2 (x y) = H 1 - ^M^l = y l ~ hx 1 ~ c f^ x iV) 
l-bfi(x,y) I -ex l-bfi(x,y) 

and hence the map 

<*■»> - ^-(^' glgglS) )- (10) 



It is straightforward to verify that 



H(x,y) = y±-^ (11) 
1 — cx 

is an invariant. 

If we denote by h the constant of the map, the reduction to one dimension yields 

X = hx^—^. (12) 
1 - fee v ' 

The other component Y is given, once X is known, by the formula Y — h(l — cX)/(l — bX). 

(|12|) is apparently nonintegrable when b = and c ^ since it is a logistic map. It is, however, integrable when 

c = or c = b since it becomes a Mobius map. Note that, if we had chosen a different fi, say x + y instead of xy, we 

obtained a nonintegrable map with the same invariant even at c = 0. We study the map 1)10(1 in detail in §4 and 5. 

For the third example we consider the case where the level set is a biquadratic curve defined by 

V, H) = ay 2 x 2 + b(y + x)yx + c(y - x) 2 + dyx + e(y + x) + f, (13) 

or equivalently 

^{x,y,H) = <j){x)y 2 +rj(x)y + p(x), (14) 

where 

4>{x) := ax 2 + bx + c, rj(x) :— bx 2 + (d — 2c)x + e, p(x) := cx 2 + ex + f. (15) 
We have not written H explicitly but assume that it is hidden in the parameters 

q = {a,b,c,d,e,f) 

which characterize the curve. 
If we fix fi(x, y) simply by 



fi(x,y)=y, (16) 



we obtain fa(x,y) from (f 1 4|) as 



hix > y) = 2^) ' (17) 



This particular choice of fi(x,y) defines a one dimensional map x — > X specified by the equation 

S(X,x;q)=0 (18) 

where 

S(X, x; q) := aX 2 x 2 + b(X + x)Xx + c(X - x) 2 + dXx + e(X + x) + f. (19) 



We shall call a map of this form a 'biquadratic map' in the following. 

The image of the map (X, Y) is determined if H is given. As an example let us assume that 

q = q' + hq", (20) 

and try to find a map which conserves H — h. Substituting 1)20)1 into 1)13(1 or 1)14)1 we see that the quantity which 
must be conserved is 

H (x v)=- ^W+jtMl+ZM (2D 
1 lU> <p"{x)y 2 +n"{x)y + p"{x)' 1 ' 

Here (</>', t/,p') and (4>", rj" , p") are defined by changing q in (0, rj,p) to q' and q" , respectively. According to our 
prescription of finding a map the image (X, Y) of (x,y) is derived by the replacement of h in (/i,/2) by 1)21)1, After 
some manipulation we obtain the map 



1 1' \v<> writ' 1 ()/.//..(-) as (.r„ -2-.''„- |. •''„.)• ' ' h nollmig 1 nit (he symiiiHric Vfvsiun ol'thc well known QRT rtjuat ionQ, 

From our construction the correspondence of this two to one map (x n , a;„+i) — > a; n +2 to the biquadratic map (|18|) is 
obvious. 



Similarly we can use the general scheme to study higher dimensional maps. For example we define 

<pi(x,y, z,i?i) = xyz-Hx, <p 2 (x, y, z, H 2 ) = (1 - x)(l - y)(l - z) - H 2 , (24) 

and fix fx by 

f , \ l-y + yz 

h{x,y,z) = x- ■ . 

1 — z + zx 

The other functions f 2 , fa are determined from (|24|) and we obtain the 3 dimensional Lotka-Volterra (3dLV) map 

l-y + yz 1-z + zx 1 - x + xy 

X = x- , Y = y- , Z — z . (25) 

1 — z + zx ' 1 — x + xy 1 — y + yz 

If we had chosen another function for fa we would have a different map but the same invariants. 
After the elimination of y and z from (|25|1 by using (|24|1 we obtain a biquadratic map of l|18|l with 

a = r + 1, b = s — 2r — 1, c = r — s, 



d = s 2 + rs + 5r - 2s + 1, e = -r(s + l), f = 0. 
In these formulae (r, s) = H 2 ) are the values of the constants of the map (|25|1 . 
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(26) 



We can carry out similar calculations for the level set 

(pi = xix 2 x 3 x 4 - Hi, ip 2 = (1 - xi)(l - x 2 )(l - x 3 )(l - Xi) - H 2 , 

ip 3 = (1 - x 2 ar 4 )(l - xix 3 ) - H 3 , (27) 

and the function 



l-x 2 + x 2 x 3 - X 2 X 3 X4, 

h{x) = X!- ■ . 

1 — Xi + X4X1 — XAX\X 2 



From (|27| we obtain the map 



1 - X2 + X 2 X 3 - X 2 X 3 X4 1 - X 3 + X 3 X 4 - X 3 X4Xx 

Xi~xi- ■ , X 2 -~x 2 - ■ , 

1 — X4 + X4X1 — XiX\X 2 1 — X\ + X\X 2 — X\X 2 X 3 

1 - £4 + X4X1 - X4X\X 2 1 - Xl + XiX 2 - X\X 2 X 3 

X 3 :=x 3 - ■ , X A :=X4- ■ , (28) 

1 - X 2 + X 2 X 3 - X 2 X 3 X4 I — X3 + X 3 X4 - X 3 X4Xi 

which is known as the discrete time Painlcvc V. The reduction of this map yields again a biquadratic map (|18J) with 
the coefficients 

a = (s + v — r + l)p + 1 — 1, 

b = (2r - s - v - 2)p - 2r - s - v + 2, 

c = (1 — r)p + r + s + v — 1, 

d = 4(1 - r)p + 2(r - l)(s + 2) + (s + v)(4 - s - v), 

e = (2r + s + v-2)p+(s + l)(v-2r-l) + (v-3)(v-l), 

f = — (r + rs + v — l)p + r + rs(r — v + 1) — (v — l) 2 , 

where (r, s, v) = (Hi, H 2 , H 3 ) and p is a solution of 

p 2 — (r — v + l)p + r = 0. 

III. DIMENSIONAL REDUCTION OF PERIODICITY CONDITIONS 

We now study, in this section, general features of periodicity conditions and present the proof of our theorem stated 
in§l. 

We consider the periodicity conditions of period n, 

X< n) = Xj, j = l,2,...,d. (29) 

When the map has p invariants, the number of independent conditions of (|29|) is d — p. Let us denote by £ = 
(£i)£2; ■■■i£,d-p) the variables which parameterize the variety V{h) of JSJ) after the reduction of p components of x. 
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The problem of finding periodic points is equivalent to finding the pth elimination ideal generated by some functions 
In satisfying 

r^(h 1 ,h 2 ,...,h p ,Ci,a,-,Cd- P )=0, a = l,2,...,d-p. (30) 

This new set of equations H3()|l plays the role of the conditions l|29l) after the dimensional reduction. 

For arbitrary values of hi, h 2 , h p , the equations (|30|1 determine a set of isolated points, hence a variety of 
dimension 0, on V(h) in general. We called in [J the periodicity conditions 'uncorrelated' in this case. There are 
possibilities that the equations r„ = impose relations on hi, h 2 , h p instead of fixing all £ Q 's. If s is the number 
of independent constraints imposed on hiS, the same number of £ a 's are left free. Accordingly the map of period n 
is extended to an s dimensional subvariety of V(h). We call the periodicity conditions 'correlated' in this case, and 
'fully correlated' when s = d — p. 

We have, so far, assumed that the values of hi, h 2 , h p have been fixed and £1, £,d- P are determined by them. 
If the periodicity conditions are fully correlated, however, none of £ a 's is constrained but only relations among the 
invariants (hi, h 2 , h p ) are imposed. To emphasize the independence from the variables £1, £2, £,d- P , in this 
particular case, we denote by j^ihi, h 2 , h p ) the d — p functions r„ (hi, h 2 , h p , £1, ^> •••) £d-p)- Note that the 
periodicity conditions are relevant only to the relations among the invariants but not to their values. This means 
that, by substituting Hi(x) for hi in Jn(h), every point on the set defined by 

««7*» = {x\ l ia) (Hi(x),H 2 (x),...,H p (x)) = 0, a= l,2,...,d- p), 

is a periodic point of period n. We called f((7n)) in Q an 'invariant variety of periodic points'. We must emphasize 
that v((j n )) is a subvariety of C d , but not of V(h). 

For the properties of v(("/ n }) we simply summarize as follows. 

1. The dimension of f((7n)) is p. 

2. Every point on v((7n)) can be an initial point of the periodic map of period n. 

3. All images of the periodic map started from a point of u((7 n )) remain on it. 

4. w((7n)) is determined by the invariants of the map alone. 

Having introduced some notions about periodicity conditions we are now going to establish the theorem stated in 
§1. To this end let us consider the periodicity conditions of period k given by 

^ a \hi,h 2 ,...,h p ) = 0, a = 1,2,..., a (31) 
T^Qii,.. ., hp, € s+ i,. £d-p) = 0, a = s + 1, ...,d-p. (32) 

They are fully correlated if s — d — p, correlated if d — p > s >0 and uncorrelated if s = 0. The fully correlated 
periodicity conditions are possible only when p > d/2, because p > s is true always. 
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We now assume that the periodicity conditions l|30|) of period n ^ k are uncorrelated. Let 
£ s -l_i(/i),£ s -l_2(/i), £d_ p (/i) be solutions of (|32[l and substitute them to i|3U|) 

T^(h 1 ,...,h p ,^,...,^,C s+ i(h),...,^- P {h))=0, a = l,...,d-p. (33) 

We consider (I33|) as a set of d~ p equations to be solved for £i, ...,£ s when the invariants are given. Since the number 
of the equations in (|33|1 is larger than the number of the unknown variables £i,...,£ s > the existence of solutions 
to l|33l) depends on the values of the invariants. From l|31|l we can decide s invariants, say hi, hi, h s , in terms 
of the rest of them. If the number p — s of the free invariants /i s +i, h p exceeds d — p — s, thus p > d/2 holds, 
we obtain solutions of l|33|) for £i,...,£ 8 , because the algebraic equations (|33|) always have solutions for the d — p 
unknown variables £1, £ s , /i s +i, h p . We now recall that the values of £1, ^2j •••> satisfying ll-i-it also satisfy 
since £s+i, £d-p are common in both and £1, £2, £s can be chosen arbitrary in 1)32)1 if s =^ 0. In other words, 
if the periodicity conditions of period n are uncorrelated, there exist points of period n which coincide with points 
of period k as long as p > d/2 and s > 1. This apparent contradiction to our assumption n ^ k shows that a set 
of uncorrelated periodicity conditions is incompatible with p > d/2 when there exists a set of correlated periodicity 
conditions. If p < d/2, on the other hand, ll-i-it are not forced to share points of other periods, hence uncorrelated 
periodicity conditions are not excluded in general. From this result our theorem presented in §1 follows immediately. 



We notice that our theorem does not exclude possibilities that all periodicity conditions are uncorrelated even 
though the map has invariants more than d/2. For example we discussed in §2 two dimensional maps © and 11UI) 
which become the logistic map after the dimensional reduction using the invariants. 



IV. SERIES OF INVARIANT VARIETY OF PERIODIC POINTS 



We have presented many examples of invariant variety of periodic points associated with some known integrable 
maps. They have been derived by solving the periodicity conditions one by one. The purpose of this section is to 
show that, for some simple cases, we can derive iteratively series of infinite number of 7„'s. 

First we study the two dimensional map 1)11) [I. After the elimination of y there is only one periodicity condition 
X( n ) = x for each period n. The variable £ that parameterizes V(h) is x in this case. The condition of period 2 is 
given by T2Q1, x) = with 

chx + bx — 1 — h 



T2(h,x) = h + 1 + chx- 



1-bx 



When c ^ 0, there are two solutions corresponding to two isolated points. When c = 0, i.e., the case of the Mobius 
map, the condition is fully correlated since 

7 2 (7i) = h + l 
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is a function of the invariant h alone. By substituting the formula for the invariant with c = we find the 
invariant variety of periodic points of period 2 as 



y{l -bx) + l = 0. 
We can continue this to the period n and obtain, in the case of c = 0, 

7nW =ln-l{h) + h n -\ 

Therefore we obtain 

h n - 1 

ln{h) = _ - , n = l,2,3, .... 
The invariant variety is (n — l)fold as it is derived by the substitution of y(l — bx) for h in 

h = e i2wk/n , fc = l,2,...,n-l. (34) 

The Mobius map, which we have just discussed, is a unique rational one dimensional map which is integrable. If wc 
consider higher dimensions we have shown in §2 that there are many rational maps which are integrable and reduce 
to the biquadratic map (|18|l of one dimension. In the rest of this section we study this map. 

The problem of finding the image Q(x) of the second iteration of the map l|18f) is equivalent to finding an elimination 
ideal generated by S(Q, X\ q) and S(X, x; q). After some manipulation we obtain an equation W 2 (Q, x) — 0, where 

W 2 {Q,x) = (Q-x) 2 S(Q 1 x ] q 2 ), (35) 

which fixes Q as a function of x. The new parameters q 2 = (ci2, b 2 , c 2 , d 2 , e 2 , f 2 ) are given by 

a 2 := (ae- cb) 2 - (ad-2ac-b 2 )(be- cd + 2c 2 ), 

b 2 := {ae-cb)(2af-be + cd-4c 2 )-(ad-2ac-b 2 )(bf-ce), 

c 2 := {af - c 2 ) 2 - (ae-bc)(bf -ce), (36) 

d 2 : = 4(0/ -c 2 ) 2 -2(ae-bc)(bf - ce) - (be - cd + 2c 2 ) 2 

-{ad-2ac-b 2 ){df - 2c/- e 2 ), 

e 2 := Ub~ce){2af -be + cd-Ac 2 )-{fd-2fc-e 2 ){ea-cb), 

h := (fb-ce) 2 -(fd-2fc-e 2 )(be-cd + 2c 2 ). 

It is not difficult to understand the meaning of the expression of the right hand side of l|35|) . The biquadratic map 
is symmetric under the exchange of the arguments, as it is manifest in l|19|) . Due to this symmetry, the backward map 
X = X( _1 )(x) solves S(X,x;q) — when the foreward map X^(x) solves it. The substitution of X = X^(x) to 
S(Q,X;q) for X yields Q — X^(x) and Q = x as its solutions, while the substitution of X = X^^fa) yields Q = x 
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and Q — X^ 2 \x). Therefore the function W2[Q,x) must be proportional to (Q — x) 2 . Moreover the elimination of 
X from S(Q, X; q) and S(X, X] q) makes W 2 {Q 1 x) quartic both in Q and x. Thence another factor in W 2 {Q^ x) must 
be symmetric and also quadratic in Q and x. The most general function satisfying these conditions is the biquadratic 
function S{Q,X]q2) itself. The parameters q 2 can be fixed by comparing both sides of (|35[1 and are given by l|36JI . 
Note that S(Q,x;q 2 ) = is solved by X {±2 ^{x). 

Now we suppose that the fcth image X^ of the map l(T%|l is found up to k — n by solving S(Q,x;qk) = 0, k = 
2, 3, n for Q. The (n + l)th image X < ~ n+1 " > of x must be obtained by solving an equation W n +i(Q, x) = which is 
derived from S(Q, X; q n ) = and S(X, x; q) = after the elimination of X. Now by the substitution of X^^x) to 
S{Q, X; q n ) = for X we obtain solutions X^ n±1 \x) and X-^ 1 \x). We know already X^^^x) as the solutions 
of S(Q, x; q n -i) = 0. Therefore W n +i(Q, x) must have S(Q, x\ q n -i) as a factor. From the same argument in the case 
of n = 2 another factor must be S(Q, x; <z„+i), hence we have 

W n+ i(Q,x) = S(Q,x]q n ^i)S(Q,x]q n+1 ). (37) 

The new coefficients q n +i will be obtained by comparing both sides of (|37|l as functions of q and q n iteratively. Using 
the notation (g/\g') n = gg' n — g'gn, the results are as follows: 



a n +i = ^—((a A c) 2 l -(a A b) n (b A c) n \ 
a n -i v ' 

b n +i = ( -((a A b) n (b A c) n - (a A c) 2 n ) + (a A c) n ( (a A e)„ + 2(b A c) n ) 

7j((a A b) n {b A e) n - (a A b) n (c A d) n + (a A d) n (b A c) r 



2 

c„+i = — ^ — [ (ce n - bf n )(ae n - bc n ) + (cb n - ea n )(fb n - ec„) 



2c„_i 



+(a/n - cc n ) + (fa n ~cc n ) , (38) 



rfn+l = j ( - /ra-lOn+1 _ On-l/n+1 - 46 n _ie n+ i - 4e„_i&„ + i + (a A /)^ + (c A d)^ 



+2 



-(a A 6) n (e A /)„ - (b A c) n (c A e) n + {a A d) n (d A f) n + 2(6 A e)„(a A /)„ 
-^3(c A e)„ - (b A f) n - {d A e) r )j ^3(& A c) n - (a A e)„ - (b A d) r )j 

{a A d) n - (a A c)„) (^(c A f) n - (d A /)„^ + 2^(6 A c)„ + (a A e)„) ((6 A /)„ + (c A e)„) J , 



7^7 (^(^W 6 ^" " ( /aC )») + (/a c )"((/a^)« + 2(e A c) n ) 
^ ^((i/\e)™(e A &) n - (f A e) n {c A d) n + (f A d) n (e A c) n 



ln+1 



7—((Iac)1 - (/ A e)„(e A c), 

fn-1 v 
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Our problem of studying the behaviour of the map (|19() has been converted to studying the behaviour of the 
parameters q n = (a n ,b n , ...,/„). These parameters determine a trajectory of iteration of the map l|18fl : 

• • • «— «— <— x — > x« — » x< 2 > 

Namely we can consider ()38|) as a map — > <Zn+i in C e . 

The periodicity conditions of period n for the map x — > X defined by (|18|l are satisfied if the parameters satisfy 

3n(9)=9- (39) 

Despite the complicated expression of the relation JSSJ, we observe a special dependence on the nth parameters q n . 
Besides c„ + i, the dependence of the [n + l)th parameters on the nth ones is always in the form (<7a<?')™ = gg' n — ff'<7n- 
They all vanish simultaneously when the periodicity conditions 1)39(1 are 'fully correlated'. In other words if there 
exists a function 7n+i(<z) such that 

q n (q) = q + j n +i(q)q n {q) (40) 

so that (|39|) is satisfied by a single condition r y n+ i(q) = 0. Here we use the same notation 7„, which we used already 
as a function of h in higher dimension. Although they arise in different contexts they turn out to be the same object 
as we will see later. 

When 140(1 holds, the equation S(Q,x;q n +i) = can be written as 

c n+1 {Q-x) 2 + 1 2 n+l {q)K n+1 {Q,x) -0. (41) 

Here 

K n+ i(Q,x) = a n+ iQ 2 x 2 + b n+ i(Q + x)Qx + d n+ iQx + e n+ i(Q + x) + f n+ i, (42) 

and d n +i, for instance, is obtained from <z n +i simply replacing (<7a<?0« ty (ffAffOn- If Q is a point of period n + 1, 
the first term of 1(41(1 vanishes. Hence the periodicity condition requires for the second term to vanish. This is 
certainly satisfied for arbitrary x if ^Jn+iiq) — 0, namely when the periodicity conditions for the parameters q n are 
fully correlated. The other possible solutions obtained by solving K n +i(x, x) — will not correspond to the points of 
period n + 1, but represent the fixed points or the points with periods which divide n + 1. 

Let us present the functions 7 n (q) explicitly in the cases of small number of n. We can show 

{a A b) 2 = (af ~eb-3c 2 + cd){2a 2 e - abd + b 3 ) 
(a A c)2 = (af — eb — 3c 2 + cd)(a 2 f + ac 2 — acd + b 2 c) 
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(&ac) 2 = (a/ -eb-3c 2 + cd)(2ace - abf -be 2 ) 



(e A /) 2 = (af-eb-3c 2 + cd)(edf~e 3 -2bf), 

from which we find 

73(g) = af -be-3c 2 + cd. (43) 

If we further calculate the cases of n = 3, 4, 

74(g) = 2ac/ - adf + b 2 f + ae 2 - 2c 3 + c 2 d - 26ce, 

75(g) = a 3 f 3 +(-cf 2 d + 2cfe 2 + fde 2 -3ebf 2 -e 4 -c 2 f 2 y 2 

+ ( - 13c 4 / + 18c 3 /d + de 3 6 + 2c/ 2 6 2 + 7dc 2 e 2 - ce 2 d 2 - 2ce 3 6 (44) 
+2c 2 feb-7fd 2 c 2 - Uc 3 e 2 + cd 3 f + /6 2 e 2 + f 2 db 2 - ebd 2 f^a 
-cd 2 b 2 f - b 3 e 3 - Ac 3 deb + cdb 2 e 2 + 13ec 4 6 - /V + 7/6 2 c 2 d 
+c i d 2 - 5c 5 d + 5c 6 - 2/6 3 ec - e 2 c 2 b 2 + eb 3 df - Ufb 2 c 3 . 

The formula l|38|l enables us to continue finding a series of~/ n (q) systematically. From the way of this construction it is 
apparent that all 7's are functions of the invariants alone if one of them is so, a fact being consistent with our theorem. 

The 3-dim Lotka-Volterra map after reduction to one dimension is a special example of the biquadratic map l|19fl . 
The parameters have been specified by (|26l) . From the general argument we obtain the result of first iteration simply 
substituting these data into Il36[l . They are given by 

<22 = (s + l) 2 s(r 2 — rs 2 — s — 3rs), 

62 = (s + l) 2 s(2r 2 s + s + 5rs - 2r 2 - r 3 - s 2 ), 

C2 = (s — r)s(r + l)(s 2 — r 2 s — 3rs — r), (45) 

d 2 = (s + l) 2 s(2rs 2 + 2s 2 - 3r 2 s - 8rs - s + r 3 s + r 4 + 5r 3 + 6r 2 - s 3 ), 

e 2 = (s + l) 2 rs(s 2 - rs + s - 2r - r 3 - 2r 2 ), 

f 2 = (s + l) 2 r 2 s(r 2 - rs + r + s 2 + s + 1). 

Notice that all parameters apart from c 2 are proportional to a common factor (s + l) 2 . This is a result which we 
cannot derive generally from H36JI since they are not factorized. On the other hand this factor is exactly the one we 
expect from our previous result [J, although this is a special behaviour of 3dLV map. 

The periodicity conditions of period 3 and higher can be derived similarly and the associated functions 73, 74 and 
75 (q) can be read off from (|43|l and l|44|l directly. We find 

73 (r, s) = r 2 + s 2 — rs + r + s + 1, 
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74(r, s) = 3rs + s + s — 3s r + r s + 6r s — r , 

7g (r, s) = r 3 s 4 - r 3 s 2 - 6r 4 s 5 + 10r 3 s 6 + 3s 5 r + s 6 + s 5 + 3r 4 s 4 - 3rV - 6r 4 s 3 
-r 6 s 3 + 3r 5 s 4 + s 4 + 21s 4 r 2 + 6s 4 r + r 3 s 7 + s 7 + 27s 5 r 2 - 3s 6 r - r 3 s 5 
+21r 2 s 6 - 10r 3 s 3 - 6rs 7 + s 8 . 

These include again precisely the conditions expected from our direct calculations in Q . 

V. TRANSITION TO INTEGRABLE REGIME 

The method of generating maps, developed in §2, enables us to study maps having invariants from a general point 
of view. For example this scheme provides us a map which becomes either integrable or nonintegrable upon reduction 
to one dimension, but shares invariants in common. Such a map will be convenient for the study of the transition 
between integrable and nonintegrable regimes. 

The map presented in (|l()(l is one such example. By adjusting properly the parameters the map can be shifted from 
nonintegrable to integrable. This is seen clearly from the expression of one dimensional map (|12|l . but not obvious 
from the expression of l|10(l written as a two dimensional map. We are going to study this particular map in this 
section in detail, since it is simple but rich enough. The purpose of this section is to explore how this transition takes 
place as c approaches zero continuously. We would like to know where the periodic points go in this limit. For this 
purpose it is sufficient to study the map l|12l) . 

Since one dimensional rational maps of degree 2 have been studied in the literature it will be convenient to discuss 
a map of the normal formal, 



after conjugating by the transformation x 



Here the parameters are related by 



1 -h' + l-h 1 ^ 
1 — h h(b — c) x 



hti = 1 + (1 - hf, 

b — c 

hence the integrable limit c = corresponds to hh' = 1 . 

Some basic data of the map (|46|l are as follows. The fixed points are at 0, oo, z p = and their multipliers are 
given, respectively, by 

2-h-h' 



ti, h, 



1 - hh' 



This map has two critical points at 



1 VI - hh' 
Z± = ~h ± ^^- 
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where the multipliers vanish. 



Repeating the map n times the image of z can be written as 

Z {n) = j z (n-l) h' + Z^ 

1 f h' + z \ ( h' + ZW \ ( hi + ZW \ ( h' + \ , , 

(49) 



h n yh- 1 + zj yh- 1 + zw J \h- 1 + z&> J yti^ + zi™- 1 ) 

If fc n -i is the degree of the rational polynomial Z^ n_1 \ we see from (|48|) that k n = 2fc n _i, hence k n = 2". Periodic 
points are obtained by solving Z^ n ' = z. When n is a prime number, the number of the periodic points of period 
n is k n — 2 = 2™ — 2. Here subtraction of 2 corresponds to two fixed points at and z p . If n is not prime but has 
divisors vx, v<z, • • • , v r they must be also subtracted. We thus find the number of periodic points of period n as 

#n = 2" - #„j - #,, 2 - • • • - #„ r - 2 

= 2™ - 2" 1 - 2^ 2^ + 2(r - 1), (50) 

which increases as fast as 

(#2, #3, #4, • ■ •) - (2, 6, 12, 30, 48, 126, 240, 504, • • •). 

In the integrable limit hh! = 1, this increase ceases owing to the cancellation of factors in the numerator and the 
denominator as seen in the second expression i|49|) . and the map becomes 

= hr n z, (51) 

It is apparent that there is no periodic point unless h satisfies ()34|) . All isolated points suddenly disappear. In order 
to uncover this trick of the transition we must know where they were right before they disappeared. After some 
manipulation we find that the right hand side of (|49|) admits the following expression 

n-2 

l[(z + h k f n ~ k ~ 2 + (hh 1 - l)P n 

z (n) _ fc=~l 

fan n—2 

Y[ (z + h k ) 2n ~ k ~ 2 + (hti - l)Q n 
fe=-i 

where P n and Q n are polynomials of z. The periodicity condition Z^ = z requires 

n-2 

(h n - 1) H (z + h k ) 2 "- k ~ 2 = (hh' - \){P n - h n Q n ) (52) 
fc=-i 

to hold. When hh 1 — 1 we obtain solutions at 

z = -i/h,-l,-h,---,-h n - 2 (53) 

if h n ^ 1. 

14 



These points are, however, certainly not periodic solution of l|5ip. They are "fossils". Suppose hh' — 1 is very small 
but not zero. Then the periodicity condition (|52|) is satisfied iff z is very close to one of Q53|). Therefore the periodic 
points were in the neighbourhood of these points right before hh' reached 1. 



The behaviour of the iteration of a map is strongly controlled by the existence of repulsive periodic points. The 
Julia set is the closure of the set of repulsive periodic points of a map, and is a source of chaotic orbits of the iteration. 
After the transition to the integrable map it is expected to disappear or become neutral points. Therefore it is quite 
useful to see how the Julia set behaves when a map is changed from nonintegrable to integrable. 

There are several equivalent ways of defining the Julia set. We adopt the following one which is convenient for our 
purpose Q. 
Definition: 

If f is a map with inverse map f~ l and Zq is a repulsive periodic point, then the Julia set J(f) of f is defined by: 



{Jf- l (zo) • (54) 



This tells us that starting from one of the repulsive points we will find all points in the Julia set simply by repeating 
the inverse map. Since this definition, applicable to one dimensional maps, requires the data of the inverse map, we 
expect to supply information complementary to our results in the previous subsection. In fact we shall prove that 
all repulsive points of the Julia set converge uniformly to the countable set of neutral points l|53|) in the integrable 
limit Q,Q. 

The inverse map of (|4f>|) is 



= -(hz-h'± y/(hz + h') 2 + 4(1 - . (55) 



Starting from a repulsive fixed point, say when h! > 1, the nth inverse map adds 2" points to J(/), and the total 
number of points in the Julia set increases as fast as 2 n+1 — 1. 

The situation, however, changes significantly when hh' = 1. In this particular case the inverse map becomes 

f-\*) = { , (56) 

Note that, when hh' = 1, — 1/h is the critical point of the map l|46(l and becomes neutral. Repeating the map n times 
starting from the maximum number of different points is n + 1, consisting of 



J n := (0, -h-\-l, -h, -h\ • • • , -h n -'j . (57) 

We have thus reproduced exactly the same series l|53f) we found before. How they could be periodic points? We can 
answer this question now, which we could not see from the observation of the forward map alone. They are periodic 
in the sense that at every step of the map i|5t)|) there is a possibility to return to the value —1/h and then continue 
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the map — 1, — h, • ■ ■ again. The Julia set in the limit of hh' — ► 1 is given by Joo and all periodic points of the map 
(I46|l must approach it. 

Now we may ask how they approach J^. This problem has been discussed in ,5] some time ago. There the following 
theorem was proved. 

All points of the Julia set of the map j^6p approach uniformly in the integrable limit. 
The proof is quoted briefly in Appendix from Q . 



VI. CONCLUDING REMARKS 

By studying periodicity conditions of a map on C with p invariants we proved our theorem presented in §1. We 
can rephrase it as 

• When p < d/2, there is no invariant variety of periodic points. 

• When p > d/2, a set of isolated periodic points and a variety of finite dimension do not exist in one map 
simultaneously. 

Because a precise notion of integrability of a map has not been known, the above statements do not provide criteria 
to distinguish an integrable map and a nonintegrable one. There are, however, some reasons to believe the following 
correspondence to hold: 

1) series of invariant variety of periodic points integrable map 

2) fractal sets of isolated periodic points nonintegrable map 

If we accept this correspondence the above statements become criteria for the integrability of a map. 

n 

We have shown, in our previous paper [lj, many examples of invariant variety of periodic points which are derived 
from known integrable maps. On the other hand a fractal set of isolated periodic points appears typically in a 
nonintegrable map and becomes a source of a chaotic trajectory of iteration of the map. The study of transition of a 
logistic map to the Mobius map, which we discussed in §5, shows an example of the above correspondence. It will 
become clear if we summarize the results of the two dimensional map l|10|) as follows: 
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periodicity conditions 


periodic points 


1-dim reduction 


c ^ 


uncorrelated 


isolated 


logistic map 


c = 


fully correlated 
j(h) = l + h+--- + h n - 1 


invariant variety 
{(x,y)\y(l - bx) = e i27rfc /"} 


Mobius map 



In the case of Hamilton dynamics with N freedom, the existence of N involutive invariants including the Hamiltonian 
guarantees integrability of the system. It will be useful to see if there is any correspondence of our theorem with that 
of Liouvillc- Arnold. If we count the momenta in the Hamilton system as independent variables, we may associate 2N 
with the dimension d of the map. The number of the invariants, which are necessary and sufficient for the Hamilton 
system to be integrable, is d/2. This is exactly the critical number of the invariants p of the map, which prohibits the 
coexistence of sets of uncorrelated and correlated periodicity conditions. 

There is, however, an important difference between a continuous time evolution and a discrete map. A set of 
uncorrelated periodicity conditions, hence a series of isolated periodic points, is possible in a map even if there are 
invariants more than N = d/2. The map (|10|l discussed above is an example. In general d — 1 invariants are not 
sufficient for a map on C d to be integrable after the reduction to one dimension. This should be contrasted with the 
fact that every first order ordinary differential equation is integrable. It is well-known that a proper continuous limit 
of the logistic map is integrable. 
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Here we quote briefly the proof of the Theorem from \^ . We can write (|55|l as 

. A(z) := hz + E(z) 
f-\z) ' 

B(z) := -h'-E{z) 



(58) 



where we defined e = hh! — 1. E(z) vanishes for small values of e. In fact we can show 



\E(z)\<R e , 



for all values of e. Here R e is given by 



u ^2 + 1 
Re = ,,, Vie 



\e\+^TT\), 



which tends to zero like y |e| as e approaches to zero. The proof of our claim follows from the facts: 
1. Using A and B defined by (|58|) the nth iteration of / _1 yields 



f- n {z) = {A U1 B U2 A U3 ■ ■ ■ B u " (z) 



V\ + Vl + ■ ■ ■ + 



(59) 



2. For any W an elememt of the form A S BW, s = 0, 1, 2, • • • , n — 1 in l|59|l lies in the neighbourhood of —h s h!: 

s-l 



\A S (BW) + h s ti 



h s E(W) + h k E{A a - k - x BW) 



k=0 



(60) 



s = 0,1,2, ...,n- 1. 

Since we assumed > 1, < 1 is satisfied as long as e is small. Hence the right hand side of i|60|) is finite 
for all s. 

3. From 1. and 2. we conclude 

r n (z )-(-h n h') < 



i 



:R f 



0, V e^O. 



i-H' 

In our present problem the Julia set is the collection of all / _n (0), hence the claim is justified. 
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